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Solar  eruptions — coronal  mass  ejections  (CMEs),  flares,  and  prominence  eruptions — have 
been  conceptualized  as  processes  consisting  of  quasi-static  build-up  and  ensuing  rapid  release 
of  stored  magnetic  energy,  usually  attributed  to  “some”  catastrophic  loss  of  equilibrium  or 
instability.  This  possible  catastrophic  transition  to  an  eruptive  state  has  been  studied  using 
idealized  two-dimensional  (2-D)  structures  such  as  linear  arcades1-3  and  flux  ropes4-'.  The 
kink  instability  is  most  often  invoked  as  the  cause  of  the  onset  of  eruptions.  The  early 
calculations  used  infinitely  long  straight  flux  ropes  as  an  idealized  model  of  large-aspect 
ratio  flux  ropes8,9,  but  extensions  to  straight  flux  ropes  of  finite  length  with  line-tied  ends 
have  been  made10,11. 

Recent  observations  and  theoretical  calculations  have  shown  that  CMEs  and  eruptive 
prominences  are  consistent  with  being  expanding  toroidal  magnetic  flux  ropes  based  on 
both  morphological  and  dynamical  properties12-14,  and  observed  CMEs  are  now  interpreted 
as  such15,16.  The  physical  processes  involved  in  the  onset  of  flux-rope  eruption  have  attracted 
renewed  interest17-19.  Kliem  and  Torok18  proposed  the  so-called  “torus  instability”  (TI)  as 
the  mechanism  that  leads  to  the  onset  and  acceleration  of  flux-rope  CMEs.  This  condition 
has  been  widely  used  to  interpret  observations  with  some  ambiguous  results19,20. 

Laboratory  experiments  were  recently  conducted  to  simulate  the  expansion  of  partial 
toroidal  pinches  with  stationary  footpoints  that  may  be  scalable  to  solar  flux  ropes21,22.  The 
former  studied  the  expansion  of  quasi-stationary  flux  ropes  through  a  background  plasma 
in  response  to  injection  of  plasma  and  poloidal  flux21,  and  the  latter  showed  that  the  flux 
ropes  can  exhibit  the  so-called  “failed”  or  “confined”  eruption22.  In  such  events  on  the  Sun, 
prominence  eruptions  stop  after  reaching,  say,  200-300  km  s_1,  becoming  quasi-stationary 
at  new  heights23,24.  In  Myers  et  al.’s  work22,  the  cessation  of  the  eruption  was  attributed  to  a 
transient  surge  in  the  tension  force  of  the  toroidal  (guide)  held.  They  compared  their  results 
to  the  TI  model  and  noted  that  the  model  lacks  a  toroidal  magnetic  held  and  its  tension  so 
that  it  cannot  model  failed  eruptions  on  the  Sun  or  in  the  experiment.  By  extension,  this  also 
raises  issues  regarding  the  TI  onset  condition  based  on  unbalancing  forces.  More  broadly, 
they  identified  this  as  an  essential  deficiency  in  the  current  theoretical  CME  models.  Indeed, 
there  is  as  yet  no  theory  to  self-consistently  model  slow  pre-eruption  evolution  leading  to 
eruption  of  CME-like  hux  ropes  under  quasi-static  build-up  of  coronal  magnetic  energy,  the 
long-standing  paradigm  of  solar  eruptions. 

This  paper  presents  the  hrst  theoretical  model — a  single  consistent  set  of  hux-rope  equa- 
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FIG.  1.  Solar  flux  rope.  The  footpoints  are  stationary  (Sf  =  const).  The  current  channel  of 
major  radius  R  and  minor  radius  a  has  toroidal  (R)  and  poloidal  (Ip)  currents,  producing  self-field 
components,  Bp  and  Bt,  respectively.  9  is  measured  from  the  footpoint  and  Z  from  Sun  center. 

tions  of  motion — to  describe  as  one  continuous  dynamical  process  the  quasi-static  flux  rope 
evolution,  catastrophic  loss  of  equilibrium  and  onset  of  eruption,  cessation  of  eruption,  and 
the  post-eruption  quasi-equilibria.  The  solutions  show  that  the  three-dimensional  (3-D)  flux 
rope  is  stable  to  kink.  Motivated  by  Ref.22,  the  forces  and  their  evolution  throughout  the 
process  are  quantified. 

Figure  1  shows  a  schematic  of  the  model  configuration  with  the  following  features.  It  is 
a  toroidal  flux  rope  embedded  in  the  corona,  with  the  two  footpoints  separated  by  S'/  = 
const.  Here,  “toroidal”  refers  to  a  section  of  a  torus.  The  flux  rope  is  defined  by  a  current 
loop  (major  radius  R  and  minor  radius  a)  and  its  magnetic  field.  The  centroid  of  the  apex 
is  at  height  Z  (from  Sun  center).  Major  radius  R{0)  is  taken  to  be  uniform  along  the  flux 
rope,  but  a(9)  varies  with  8\  at  the  apex,  a{0  =  7r)  =  aa,  and  at  the  footpoints,  a{0f )  =  a/. 
If  unsubscripted,  a  =  aa  is  understood. 

The  current  channel  carries  toroidal  and  poloidal  currents,  R  and  Ip,  producing  Bp  and  Bt , 
respectively,  via  Ampere’s  law  The  current  is  localized  within  the  channel,  i.e.,  J(r  >  a)  —  0, 
where  r  is  the  minor  radial  coordinate.  Thus,  R  =  2i r  f"  Jt(r)rdr  and  Ip  =  2itR  JQa  Jp(r)dr, 
with  Bpa  =  Bp[r  =  a)  =  2 It/ca.  Outside  the  current  channel  (r  >  a),  Bp(r )  ~  Bpaa/r.  By 
fiat,  the  flux  rope  is  taken  to  be  limited  to  r  <  2a,  beyond  which  Bp{r )  ~  Bpa/2  may  be 
too  weak  to  be  distinguished  from  the  ambient  field.  Representative  magnetic  field  “lines” 
are  drawn.  The  toroidal  flux  is  <f>t  =  27t  J"  Bt(r)rdr  =  tt Bta2,  and  —  const  is  assumed, 
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where  Bt  is  the  minor  radial  average. 

The  flux  rope  has  average  internal  pressure  p,  embedded  in  the  corona  of  pressure  pc(Z ) 
and  magnetic  field  BC(Z).  Here,  Bc  is  perpendicular  to  the  plane  of  the  flux  rope.  These 
system  parameters  have  been  defined  elsewhere2'.  There  can  be  an  external  toroidal  com¬ 
ponent  Bct,  which  will  not  be  treated  because  of  space  limitation  and  the  need  to  introduce 
a  spatial  model  Bct(Z). 

The  ideal  MHD  force  density  is  given  by  f  =  (l/c)J  x  B  —  Vp  +  pV</>9,  p  is  mass  density, 
and  (j)g  is  the  gravitational  potential.  The  volume-integrated  major  radial  force  per  unit 
length  at  the  apex  is  given  by25,26 

m1F  =  ?Mif!c(t)  +  Fd  +  Fs’  (1) 


where  M  =  na2p  is  the  mass  per  unit  length,  and 


B2t 
2  B2a 


+ 


2  R  Br 


a  B, 


pa 


1  +  -. 
2 


(2) 


Here,  f3p  =  8ir(p  —  pc)  /  B2a  is  pressure  gradient  (dp /dr)  integrated  over  a.  The  minor  radius 
a(t )  is  governed  by 


d2a  I? 

M - =  —  , 

dt 2  c2a  V  B‘i 


B2 


—  i  +  A 


pa 


(3) 


In  equation  (1),  $p(f)  =  cL(t)It(t )  is  the  poloidal  flux  of  the  flux  rope,  so  that  <h2/c4L2  = 
/2/c2.  Neglecting  Fd  and  Fg  for  now,  the  equation  can  be  written  as 


d2Z_R 
~dF  ~  ^ 


(4) 


R 


where  rR  =  R/Vav  with  Vav  =  Bpa/ (ArK~p)1^2 .  Thus,  the  acceleration  timescale  is  the  Alfven 
transit  time  tr.  For  the  main  acceleration  phase,  rR  depends  on  the  initial  equilibrium 
values  of  p0  and  Bpa0. 

The  toroidal  geometry  enters  the  equations  via  the  constraint 

Z2(t)  +  S2/A 


R(t)  = 


2  Z(t) 


(5) 


where  Z  =  (Z  —  f?Q).  That  is,  at  any  time,  R(t)  is  determined  by  Sf  =  const  and  apex 
height  above  the  photosphere  Z(t ).  The  self-inductance  L(t)  relates  the  current  It  to  the 
poloidal  energy  Up  =  (1/2 )LI2  and  is  given  by2' 


L(t)  = 


AttQR 


lnf^+mf-) 

V  af  J  \  aa  J 


2+  -  +  A 
2 


(6) 
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where  27 rQ(t)R(t)  is  the  length  of  the  evolving  flux  rope  and  0(f)  is  the  fraction  of  a  complete 
circle.  The  additive  quantity  A (f)  is  a  geometrical  factor  and  is  determined  by  how  a{6) 
increases  from  aj  to  aa,  with  A (f)  =  0  for  a{6)  =  const  or  a{6)  ~  exp[(0  —  Of)/s],  where 
s  =  const.  We  will  use  exponential  dependence  so  that  A (f)  =  0.  For  CME-likc  flux  ropes 
with  linear  a(9),  |A|  ~  0. 1-0.3. 

The  drag  force  per  unit  length  is  modeled  by 

Fd  =  4cdPcaAV\AV\,  (7) 

where  AV  =  Vc  —  (V  +  2 w),  V  =  dZ/dt  is  the  speed  of  the  centroid  of  the  apex  at  Z(t), 
w(t )  =  da(t)/dt,  and  Cd  is  the  drag  coefficient.  The  term  ( V  +  2 w)  is  the  speed  of  the 
outermost  magnetic  surface  of  the  flux  rope  (r  =  2a)  and  arises  from  the  definition  of  the 
leading  edge  (LE),  ZLE  =  Z  +  2a  (Figure  1),  so  that  Vle  =  d(Z  +  2a) /dt  =  V  +  2 w.  The 
values  of  cd  suitable  for  coronal  and  solar  wind  applications  have  been  found  to  be  in  the 
range  of  0.1-3  in  an  MHD  simulation  study28,  and  cd  =  1  is  used  in  the  calculations  here. 

The  terms  in  fE  were  originally  derived  by  Shafranov29  for  axisymmetric  tokamak  equi¬ 
libria  in  vacuum.  Equations  (1)  and  (3)  were  applied  to  the  dynamics  of  non-axisymmetric 
solar  flux  ropes  with  Sf  =  const  and  constraint  (5),  including  drag  coupling  to  the  ambient 
corona26.  The  hoop  force  for  arbitrary  a/R  <C  1  was  later  rigorously  derived30.  The  hoop 
force  is  now  familiar  in  the  solar  context 18,22,25,28’31. 

The  initial  coronal  loop  is  defined  by  footpoint  separation  Sf,  apex  height  Zq,  and  aspect 
ratio  Ro/ao,  with  average  internal  pressure  po-  Major  radius  Rq  is  calculated  from  Sf  and  Z0. 
The  structure  is  in  equilibrium  with  coronal  pressure  pc o  =  pc(Z0 )  and  overlying  magnetic 
field  Bc0  =  BC(Z0).  Given  the  initial  geometry  and  coronal  parameters  pc0  and  Bc0,  the 
equilibrium  force  balance  conditions,  d2Z/dt2  =  0  and  d2a/dt2  =  0,  are  imposed.  These 
conditions  determine  the  initial  equilibrium  field  Bpa o  =  Bp{r  =  a\t  =  0) 

Bpao  =  2  tAo  +  fdpo  ~  4:7rrn,Rog0(no  -  nc0)]  \Bc0\ ,  (8) 

so  that  Bpa0  oc  (R0/a0)\Bc0\  >  \Bc0\.  Here,  A0  =  ln(8i?0/ao)  —  (3/2)  +  (£j/2).  Equation  (3) 
yields  the  toroidal  field  B20  =  (1  —  /3po)B2a0.  Thus,  Bp  and  Bt  both  critically  depend  on  the 
toroidal  geometry  ( R/a ).  With  Bc  ^  0,  it  is  possible  to  choose  a  force-free  initial  equilibrium. 
Thus,  we  will  use  (3po  =  0  in  which  case  Bt o  =  Bpa 0.  The  structure,  however,  is  not  assumed 
to  evolve  force-free,  i.e.,  /3p(t )  ^  0  for  t  >  0. 
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FIG.  2.  Quasi-static  evolution,  (a)  Profiles  of  d$p(t)/dt.  The  reference  profile  (solid  curve)  is 
labeled  “1.”  Respectively,  profiles  “1,”  “2,”  and  “3”  have  ( d&p/dt)max  =  (5.8, 7.0, 4.2)  x  1015 
Mx  s-1].  (b)  Height  of  the  apex  Z(t)  Zo  =  I.IIRq  is  indicated,  (c)  Speed  of  the  apex,  V  =  dZ/dt. 
Solutions  “1,”  “2,”  and  “3”  for  the  three  d$p(t)/dt  profiles  in  panel  (a). 

The  discussion  below  will  treat  a  CME-like  flux  rope  to  illustrate  the  theory  in  real  units 
so  that  the  theoretical  results  may  be  directly  related  to  observations.  The  initial  flux  rope 
is  specified  by  S'/  =  1  x  105  km,  Z0  =  8.0  x  104  km,  and  Ro/clq  =  2  (a0  =  2.8  x  104  km), 
giving  Rq  =  5.6  x  104  km.  The  coronal  parameters  at  Z  =  Zq  are  Bc o  =  —4.0  G  and 
pc o  =  0.32  dyn  cm-2.  With  T0  =  Tc0  =  2  x  106  K,  we  have  h(0)  =  5.8  x  108  cm-3.  These 
values  yield  Bpa( 0)  =  Bt0  =  8.64  G,  so  that  rR  =  Rq/Vapo  =  72  sec  with  Vap o  =  783  km  s-1. 
The  total  mass  of  this  flux  rope  is  MT  =  1.1  x  1015  g.  The  initial  poloidal  flux  is  <3>p(0)  = 
7.4  x  1020  Mx  (cgs)  and  =  2.1  x  1020  Mx. 

We  now  slowly  increase  4>p,  i.e.,  the  magnetic  field  “twist,”  by  specifying  d$p(t)/dt  as 
shown  in  Figure  2a.  The  maximum  rate  is  (d&p/dt)max  =  5.8  x  1015  Mx  s^1,  which  is 
sustained  for  approximately  40  hours,  with  <1 )p{t  =  60)  ~  2.34>po-  This  is  designated  as  the 
reference  profile  (solid  curve  “1”).  The  injection  rate  is  “slow”  on  the  acceleration  timescale 
tr,  satisfying  ( d&p/dt)max  (&p0/tr)  ~  1019  Mx  s^1.  As  $p(t)  increases,  the  terms  in  fR 
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FIG.  3.  Potential  barrier  due  to  BC(Z)  (solid)  and  critical  height  Zcrt  =  1.18I?©  of  loss  of  equilib¬ 
rium.  Apex  velocity  V(Z)  (dashed). 

become  unbalanced  (fn  >  0),  and  the  flux  rope  rises.  Panels  2b  and  2c  show  the  evolution 
of  the  flux  rope.  The  solution  shows  that  the  apex  of  the  flux  rope  Z  slowly  rises  until  about 
t  =  36  hrs  when  it  starts  to  erupt.  This  eruption,  however,  “fails”  to  sustain  itself  after 
reaching  the  peak  speed  of  V  ~  200  km  s-1.  The  apex  velocity  V ( t )  for  this  solution  is  shown 
in  panel  c  (curve  “1”).  The  value  of  <3>p  remains  virtually  unchanged  during  the  eruption. 
Post-eruption,  the  flux  rope  reverts  to  quasi-equilibrium  evolution  at  greater  heights  and 
tends  to  a  new  equilibrium  as  d$p/dt  — >■  0.  The  full-width-at-half-maximum  (FWHM) 
duration  of  the  velocity  peak  is  At  ~  30  min. 

Another  solution  was  obtained  for  the  same  initial  flux  rope  subject  to  a  faster  rate  of 
flux  injection  (“2”,  panel  a):  ( d&p/dtd)max  =  7.0  x  1015  Mx  s-1.  The  pre-eruption  rise  is 
slightly  faster,  reaching  the  same  critical  point  earlier.  The  apex  speed  is  shown  in  Figure  2c 
(dashed,  “2”).  Except  for  the  shift  in  time,  the  profile  of  V(t )  is  indistinguishable  from  that 
of  solution  1.  The  eruption  is  insensitive  to  the  value  of  d&p/dt  provided  ( d&p/dt)max  <C 
&po/tr.  This  is  because  tr  does  not  depend  on  d$p/dt.  Profile  “3”  (dotted)  has  a  slower 
rate  (d$p/dtd)max  =  4.2  x  1015  Mx  s-1.  The  apex  velocity  profile,  labeled  “3”  (dotted)  in 
Figure  2c  is  again  unchanged  except  that  Zcrt  is  reached  later.  The  same  profile  of  failed 
eruption  is  obtained  as  the  rate  (d$p/dtd)max  is  made  smaller  (not  shown). 

Figure  3  shows  the  magnitude  of  the  external  magnetic  held  \BC(Z)\  (solid  curve),  which 
produces  a  downward  force  ( ItBc )  so  that  the  peak  in  |  BC(Z)  |  represents  an  effective  potential 
barrier.  The  initial  apex  position  of  the  flux  rope,  Z$  l.lli?Q,  is  indicated  by  the  open 
circle.  This  shows  that  the  flux  rope,  which  is  stable  to  major  radial  expansion,  must 
overcome  the  barrier.  As  <f>p  increases  slowly,  the  flux  rope  rises  (Figure  2b).  The  solid 
circle  indicates  Zcrtj  height  of  the  “onset”  of  the  eruption,  which  is  arbitrarily  chosen  to  be 
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FIG.  4.  Dynamics  during  the  failed  eruption,  (a)  V(t)  for  solutions  “1”  (solid)  and  “2”  (dashed). 
Curve  “1”  magnified  by  100  (dotted).  Diamond:  V  =  1  km  s  .  (2)  Forces.  F^/:  hoop  force.  FBc: 
ItBc.  Ftens:  tension,  the  B/  /  B/a  term  in  equation  (2).  The  net  force  Ftot  is  magnified  by  50. 

when  V  =  5  km  s'1  is  reached.  For  solution  1,  Zcrt  =  1.20 i?0,  which  occurs  at  Tcrt  =  36.1  hr. 
To  show  the  relationship  of  V  to  BC(Z),  the  apex  velocity  is  shown  as  a  function  of  Z  (dashed, 
right  vertical  axis).  The  onset  height  Zcrt)  well  past  the  peak  of  \BC\,  is  the  same  for  all 
three  solutions  shown  in  Figure  2c. 

The  apex  velocity  V (t)  and  the  forces  on  the  flux  rope  are  shown  on  an  expanded  timescale 
in  Figure  4.  The  curves  marked  “1”  (solid)  and  “2”  (dashed)  are  the  respective  curves  in 
Figure  2c.  Tcrt  for  solution  1  is  marked.  Figure  4a  also  shows  velocity  curve  “1”  multiplied 
by  100  (dotted).  The  open  diamond  designates  the  time  before  which  V  <  1  km  s-1.  This 
is  much  earlier  than  t  =  Tcrt.  Figure  4b  shows  the  various  force  terms  in  equation  (1):  the 
hoop  force  Fhj  arises  from  JtBp,  tension  Ftens  oc  —Bt/Bpa  from  the  self  held  Bt,  FBc(<  0) 
from  ItBc  term,  and  Fp  oc  /3P  from  dp/ dr  ^  0.  The  net  major  radial  force  is  denoted  by  Ftot. 
The  initial  equilibrium  (Ftot  =  0)  is  established  by  Fhf  =  ~(FBc  +  Ftens ),  with  Ftens  playing 
a  relatively  minor  role.  Because  (3po  =  0  is  chosen,  Fp0  —  0  at  t  —  0.  Except  during  the 
eruption,  Fd  is  negligible  because  AV  ~  0  in  equation  (7).  Gravity  is  negligible. 

Figure  5  shows  forces  on  an  expanded  vertical  scale.  The  net  force  Ftot  driving  the  failed 
eruption  is  a  relatively  small  imbalance  between  FBc  and  Fhf  (Figure  4b),  but  the  post- 


FIG.  5.  Evolution  of  forces.  All  forces  are  shown  in  real  units.  Ftot  includes  all  terms  in  equation  (1), 
coupled  to  equation  (3).  \Fp\  increases  during  the  eruption,  establishing  a  new  non- force- free 
equilibrium,  Fp  ^  0.  Fbc  rapidly  vanishes  before  t  =  36.6  and,  for  clarity,  is  not  shown. 

eruption  equilibria  are  established  by  Fhf  =  —  (Fp  +  Ftens)  with  comparable  Fp  and  Ftens 
and  FBc  cs  0.  (For  clarity,  FBc  is  not  shown;  it  becomes  negligible — \FBc\  <C  \Ftens  +  Fp | — by 
t  ~  36.6  hrs.)  Although  the  initial  flux  rope  is  force-free  (FJxB  =  Fp  =  0),  it  becomes 
strongly  non-force-free  (Fp  ~  Ftens).  Thus,  it  is  essential  to  include  the  f3p  term  in  the 
equation  of  motion.  Note  that  Fp0  =  0  at  t  =  0  is  possible  because  Bc  ^  0.  If  Bc  =  0,  then 
Fp  7^  0  is  necessary  for  equilibrium32. 

The  expansion  prior  to  and  after  the  failed  eruption  is  quasi-static:  it  can  be  made 
arbitrarily  slow  by  taking  the  limit  ( dQp/dt)max  — >■  0  and  increasing  the  duration  of  flux 
injection.  The  three  solutions  in  Figure  2c  illustrates  this  limiting  process.  Given  BC(Z)  and 
pc{Z),  the  velocity  and  timescale  rB  are  independent  of  how  the  limit  is  taken  (Figure  4a), 
being  determined  by  the  intrinsic  flux  rope  properties  such  as  the  geometry  and  self  magnetic 
field.  The  transition  to  eruption  is  catastrophic  in  that  Tr  is  finite. 

Numerical  simulations  of  the  evolution  of  arcades  under  “slow”  driving  of  footpoints 
have  been  carried  out33-35.  In  such  studies,  the  lower  boundaries  are  placed  at  the  base  of 
the  corona,  and  for  numerical  tractability,  the  footpoints  are  driven  much  faster  than  the 
photospheric  speeds  of  observed  magnetic  features  (e.g.,  up  to  tens  of  kilometers  per  second 
in  the  above  simulations).  To  the  extent  that  large  coronal  structures  and  photospheric 
magnetogram  features  (e.g.,  active  regions)  nearly  co-rotate,  such  high  speeds  on  the  scale 
of  eruptive  arcades  would  require  comparably  fast  speeds  the  photosphere.  These  simulations 


9 


FIG.  6.  Instantaneous  value  of  fl2(t)  =  u2/Vjpk2  (solid  curve).  Kadomtsev  condition  for  kink  sta¬ 
bility.:  Q2  >  0  means  stable  oscillatory  motion.  Characteristic  magnetic  pitch  F(t )  =  Bpa(t)  /  Bt(t) 
(dashed  curve).  The  values  are  shown  on  the  right  vertical  axis. 


show  coronal  response  on  the  Alfven  transit  time  similar  to  Tr.  In  these  simulations,  the 
onset  of  eruption  depends  on  various  “tether  cutting”  processes  via  reconnection,  which 
dynamically  reduces  the  effective  potential  barrier. 

To  examine  the  kink  stability  of  the  flux  rope  with  no  guide  field,  we  apply  the  Kadomtsev 
condition  based  on  an  energy  principle,  which  yields  the  eigenfrequency36 


ca 


2 


1  -r2(t)  In 


(9) 


where  is  T(t)  =  Bpa(t)  /  Bt(t)  is  the  characteristic  magnetic  pitch  and  k  the  toroidal  wavenum¬ 
ber.  The  flux  rope  is  stable  if  uj2  >  0.  For  the  longest  wavelength  between  the  stationary 
footpoints,  we  have  1/ka  =  QR/a.  Figure  6  shows  Q2(t)  =  c o2  jV\pk2  >  0  (solid).  This 
shows  that  the  flux  rope  is  stable  to  the  kink  (fl2  >  0)  at  t  =  0  and  throughout  the  quasi¬ 
equilibrium  expansion.  The  flux  rope  evolves  to  a  more  stable  configuration  (more  positive 
fl2)  during  the  failed  eruption.  There  is  no  external  toroidal  guide  field  ( B  ct  —  0). 

Also  plotted  in  Figure  6  is  magnetic  pitch  T(f)  (dashed);  it  decreases  during  the  eruption 
but  is  otherwise  virtually  constant.  This  is  equivalent  to  tension,  Bt{t)  /  Bp(t)  =  r_1(f)  ~ 
QR/a ,  increasing  relative  to  the  hoop  force.  Both  render  the  flux  rope  even  more  kink  stable. 
These  two  effects  are  coupled,  due  to  a  number  of  dynamical  effects:  R(t)  and  therefore  L(t) 
increase  so  that  It(t)  oc  1  / L  ~  1  / R{t)  decreases;  a{t)  expands  rapidly,  causing  Bp  oc  It/ a  and 
Bt{t )  oc  <Ib/a2(f)  to  decrease;  Fhf  does  work  against  Ftens  during  the  expansion,  equation  (1), 
so  that  Bt  gains  energy  from  Bp.  The  expansion  in  a{t )  causes  p{t)  to  decrease  independently 
of  pc{Z),  and  Fp  oc  j3p  becomes  nonzero  (non-force-free).  These  dynamical  quantities  evolve 
self-consistently  according  to  equations  (1)  and  (3),  which  determine  the  inductance  L(t) 
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and  therefore  the  evolution  of  the  magnetic  field.  This  is  a  fully  3-D  ideal  MHD  process. 

The  above  discussion  shows  that  the  minor  radial  dynamics  a(t)  are  fundamental  to  the 
evolution  of  major  radial  forces.  Figure  7  gives  the  dynamics  of  a(t)  and  R(t)/a(t),  where  a(t) 
and  R{r )  are  determined  by  coupled  equations  (l)-(6).  It  shows  that  a(t)  (solid)  rapidly 
expands  during  the  failed  eruption  while  R(t)/a(t )  (dashed)  undergoes  a  slight  decrease 
when  V(t)  peaks:  a(t)  expands  faster  than  R(t)  during  this  time.  Past  this  point,  a(t)  and 
R(t)/a(t )  both  increase,  showing  that  major  radial  expansion  dominates  again.  Throughout 
the  process,  R(t)/a(t )  7^  const.  This  is  in  contrast  to  the  assumption  of  R(t)/a(t )  =  const 
in  the  TI  model18,  which  is  not  made  in  the  calculations  of  Ref.31.  The  nontrivial  variation 
in  R(t)/a(t )  shows  that  it  is  essential  to  solve  the  coupled  equations  (1)  and  (3). 

The  eruption  described  above  is  consistent  with  failed  eruptions  on  the  Sun23,24.  The 
velocity  profile  in  Figure  4a  resembles  the  observed  profile  reported  in  Figure  2e  of  Ref.24, 
including  a  damped  oscillation.  The  FWHM  duration  At  of  the  velocity  peak  is  several 
minutes23  to  a  few  tens  of  minutes24.  The  failed  eruption  was  interpreted  as  “probably  due 
to  the  strong  strapping  field.”  The  strength  of  this  held  is  not  measured  but  should  be 
weaker  at  the  new  height.  For  the  model  flux  rope,  the  strapping  held  Bc  vanishes,  and  it 
is  Bt  and  dp/ dr  that  provide  the  equilibrium. 

For  the  hux  ropes  in  the  experiment22,  Myers  et  al.  analyzed  the  forces  in  terms  of  inte¬ 
grated  quantities  as  in  equations  (1)  and  (3).  In  one  regime,  the  hux  ropes  stabilized  against 
the  kink  by  a  sufficiently  strong  toroidal  guide  held  exhibited  failed  eruption,  returning  to 
nearly  the  original  height  after  the  eruption.  They  observed  a  transient  surge  in  the  toroidal 
guide  held  and  its  tension  force  during  the  eruption,  which  was  identified  as  the  cause  of  the 
failed  eruption.  Based  on  their  analysis  and  the  lack  of  toroidal  helds  in  the  TI  equation, 
Myers  et  al.  concluded  that  the  TI  model — by  extension  other  models  with  no  toroidal 
helds — cannot  explain  the  laboratory  or  observed  failed  eruptions,  or  the  onset  condition. 
The  TI  equation  is  a  special  case  of  equation  (1)  with  (3p  =  0  and  Bt  =  0.  In  this  limit, 
both  Ftens  and  Fp  are  absent,  the  two  forces  that  determine  the  failed  eruption  and  subse¬ 
quent  equilibria.  This  is  also  rehected  in  the  fact  TI  solutions  do  not  exhibit  any  oscillations 
that  are  often  observed21.  The  forces  and  the  length  Sf  =  const  in  the  experiment  better 
correspond  to  those  in  equations  (l)-(5)  than  the  TI  equation. 

It  is  noted  that  no  external  guide  held  ( Bct )  or  strapping  held  ( Bc )  is  required  to  confine 
the  eruption:  the  forces  involved  in  the  failed  eruption  and  the  subsequent  non-force-free 
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FIG.  7.  Dynamics  of  minor  radius  a(t)/ao  and  aspect  ratio  R(t)/a(t). 

equilibrium  discussed  herein  are  due  to  self  fields  and  pressure  gradient. 

The  results  of  the  present  work  are  consistent  with  the  argument  of  Myers  et  al.  that 
toroidal  tension  force  is  the  cause  of  failed  eruptions22.  There  are,  however,  significant 
differences.  The  model  system  has  no  conducting  vessel  nor  an  external  magnetic  held.  The 
self  Bt  and  its  tension  exhibit  no  surge  during  the  eruption,  and  the  equilibrium  flux  rope 
organized  by  self  fields  and  ambient  pressure  is  kink  stable.  Indeed,  pressure  force,  Fp  oc  (3p — 
not  considered  in  TI  and  interpretations  of  observations  and  the  experiment  results — is 
comparable  to  Ftens.  In  the  experiment,  a  conversion  of  <f>p  to  was  inferred.  This  implies 
the  presence  of  non-ideal  or  inductive  effects.  In  the  theory,  there  is  no  flux  conversion 
because  =  const,  but  Bp  does  work  on  Bt.  The  surge  detected  in  the  experiment  is 
attributed  to  the  interaction  between  the  expanding  flux  rope  and  the  applied  guide  held. 
This  interaction  may  involve  the  conducting  wall  (e.g.,  image  currents)  and  external  circuit. 
Assessment  of  such  laboratory  factors  is  identified  as  an  open  question  for  further  work. 

There  is  an  ambient  medium  (pc  ^  0)  so  that  /3P  <  0  in  the  post-eruption  equilibria.  The 
dynamics  shown  in  Figure  4  are  critically  determined  by  the  self  helds  Bt(t),  Bp(t),  and 
pressure  gradient  dp(t)/dr ,  with  the  drag  force  being  significant  only  during  the  eruption. 
The  model  hux  rope  reaches  a  new  equilibrium  height  while  the  hux  ropes  in  the  experiment 
return  to  the  initial  positions.  In  observed  failed  eruptions,  prominences  typically  do  not 
return  to  the  initial  heights. 

The  TI  model  defines  the  onset  condition  for  an  axisymmetric  torus  with  no  footpoints 
using  the  so-called  decay  index  n  =  —\(Z/Bc)dBc/dZ\crt.  In  the  experiment,  the  failed 
eruptions  occupy  the  n-q  parameter  space  region,  0.8  <  n  <  4  (TI  unstable)  and  0.8  <  q  < 
1.8  (kink  stable),  where  q  =  ( a/QR)Bt/Bpa  is  the  safety  factor.  Because  of  the  different  form 
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of  BC(Z)  and  force  terms  in  the  TI  equation,  there  are  no  directly  comparable  quantities 
here.  Nevertheless,  formally,  the  above  quantities  are  q  ~  0.5  and  n  ~  2  for  the  present 
flux  rope.  Here,  the  characteristic  value  of  q  is  based  on  Bt  because  there  is  no  external 
guide  held.  The  TI  prediction  is  ncr  ~  0.8.  The  difference  arises  from  the  toroidal  magnetic 
tension  force  Ftens  due  to  the  self  held  Bt.  The  onset  condition  given  by  the  more  general 
treatment31,37,  having  Fp  ^  0,  Ftens  ^  0,  and  Sf  =  const ,  are  applicable  to  the  model  and 
solar  hux  ropes. 

An  observed  failed  filament  eruption23  was  simulated  using  an  MHD  model38.  The  initial 
structure  is  a  kink-unstable  hux  rope  with  hxed  footpoints.  There  is  an  applied  toroidal  held 
specihed  by  a  hetitious  subphotospheric  line  current.  The  structure  is  assumed  to  be  force- 
free.  The  simulation,  with  increased  diffusion  to  prevent  numerical  instabilities,  shows  that 
the  kink  instability  grows  and  then  saturates,  which  is  attributed  to  (numerical)  reconnection 
above  the  hux  rope.  No  evolution  to  the  unstable  initial  state  is  given.  The  failed  eruption 
here  is  non- force- free  and  kink  stable,  and  requires  no  reconnection.  Nevertheless,  the  major 
radial  forces  (without  the  external  held)  should  be  describable  by  equations  (1)  and  (3). 

There  are  several  straight  2-D  models  of  catastrophic  transition0  1 .  They  typically  use 
force-free  (Fp  =  0)  cylindrical  hux  ropes  with  no  internal  toroidal  held,  in  which  the  upward 
force  arises  from  the  repulsion  from  photospheric  image  currents.  In  3-D  hux  ropes,  the 
force  from  the  photospheric  return  current  is  insignificant  relative  to  the  hoop  force  unless 
Z/ Sf  -C  l30.  There  is  no  2-D  counterpart  of  the  major  radial  forces. 

The  paper  has  presented  the  hrst  theoretical  calculation  of  quasi-static  evolution,  catas¬ 
trophic  loss  of  equilibrium,  failed  eruption,  and  re-establishment  of  new  toroidal  equilibrium. 
The  driver  of  the  evolution  is  a  quasi-static  increase  in  the  poloidal  magnetic  energy  (oc  <L2). 
Physically,  this  may  arise  from  slow  differential  twist  in  the  footpoints  or  changes  in  the  cur¬ 
rent  below  the  photosphere.  The  key  findings  are:  (1)  an  initial  force-free  flux  rope  evolves 
into  a  non-force-free  flux  rope  maintaining  kink  stability  and  (2)  the  onset  of  eruption,  cessa¬ 
tion  of  the  eruption,  and  the  post-eruption  equilibria  are  critically  determined  by  the  toroidal 
self  held  and  pressure  gradient,  i.e.,  by  Fhf  —  ( Ftens  +  Fp),  which  have  been  neglected  in  pre¬ 
vious  theoretical  models.  The  rapid  decrease  in  Fhf  is  governed  by  L(t )  oc  R(t)  \n(R / a f) . 
The  model  dynamics  are  idealized  but  are  fully  3-D,  with  equations  (l)-(8)  coupling  the 
features  in  Figure  1. 

The  work  was  supported  by  the  Naval  Research  Laboratory  Base  Research  Program. 


13 


REFERENCES 


1Low,  B.  C.,  Astrophys.  J.,  212,  234  (1977). 

2Jockers,  K.,  Solar.  Phys.,  56,  37  (1978). 

3Finn,  J.  M.,  and  Chen,  J.,  Astrophys.  J.,  349,  345  (1990). 

4van  Tend,  W.,  and  Kuperus,  M.,  Solar  Phys.,  59,  115  (1978). 

5Forbes,  T.  G.,  J.  Geophys.  Res.,  95.  11919  (1990). 

6Forbes,  T.  G.,  and  Isenberg,  P.  A.,  Astrophys.  J.,  373,  294  (1991). 

7Isenberg,  P.  A.,  Forbes,  T.  G.,  and  Demoulin,  P.,  Astrophys.  J.,  417,  368  (1993). 
8Sakurai,  T.,  Publ.  Astron.  Soc.  Jpn,  28,  177  (1976). 

9Hood,  A.  W.,  and  Priest,  E.  R.,  Solar  Phys.,  64,  303  (1979). 

10Velli,  M.,  Einaudi,  G.,  and  Hood,  A.  W.,  Astrophys.  J.,  350,  428  (1990). 

11Ryutov,  D.  D.,  Cohen,  R.  H.,  and  Pearlstein,  L.  D.,  Phys.  Plasmas,  11,  4740  (2004). 
12Chen,  J.,  Howard,  R.  A.,  Brueckner,  G.  E.,  Santoro,  R.,  Krall,  J.,  Paswaters,  S.  E.,  St. 
Cyr,  O.  C.,  Schwenn,  R,,  Larny,  P.,  and  Simnett,  G.  M.,  Astrophys.  J.  Lett.,  490,  L191, 
(1997). 

13Wood,  B.  E.,  Karovska,  M.,  Chen,  J.,  Brueckner,  G.  E.,  Cook,  J.  W.,  and  Howard,  R.  A., 
Astrophys.  J.,  512,  484  (1999). 

14Krall,  J.,  Chen,  J.,  Duffin,  R.  T.,  Howard,  R.  A.,  and  Thompson,  B.  J.,  Astrophys.  J., 
562,  1045  (2001). 

15Thernisien,  A.  F.  R.,  Howard,  R.  A.,  and  Vourlidas,  A.,  Astrophys.  J.,  652,  763  (2006). 
16Krall,  J.,  Astrophys.  J.,  657,  559  (2007). 

17T5rok,  T.,  Kliem,  B.,  and  Titov,  V.  S.,  Astron.  Astrophys.,  413,  L27  (2004). 

18Kliem,  B.,  and  Tdrok,  T.,  Phys.  Rev.  Lett.,  96,  255002  (2006). 

19Demoulin,  P.,  Vourlidas,  A.,  Pick,  M.,  and  Bouteille,  A.,  Astrophys.  J.,  750,  147  (2012). 

Erratum:  Astrophys.  J.,  754,  156  (2012). 

20Liu,  Y.,  Astrophys.  J.  Lett.,  679,  L151  (2008). 

21Tripathi,  S.  K.  P.,  and  Gekelman,  W.,  Sol.  Phys.,  286,  479  (2013). 

22Myers,  C.  E.,  Yarnada,  M.,  Ji,  H.,  Yoo,  J.,  Fox,  W.,  Jara- Almonte,  J.,  Savelieva,  A.,  and 
DeLuca,  E.  E.,  Nature,  528,  526  (2015). 

23 Ji,  H.,  Wang,  H.,  Schmahl,  E.  J.,  Moon,  Y.-J.,  and  Jiang,  Y.,  Astrophys.  J.,  595,  L135 
(2003). 


14 


24Cheng,  X.,  Hao,  Q.,  Ding,  M.  D.,  Liu,  K.,  Chen,  P.  F.,  Fang,  C.,  and  Liu,  Y.  D.,  Astrophys. 
J.,  809,  46  (2015). 

25Chen,  J.,  J.  Geophys.  Res.,  101,  27499  (1996). 

26Chen,  J.,  Astrophys.  J.,  338,  453  (1989). 

27Chen,  J.,  and  Garren,  D.  A.,  Geophys.  Res.  Lett.,  20,  2319  (1993). 

28Cargill,  P.  J.,  Chen,  J.,  Spicer,  D.  S.,  and  Zalesak,  S.  T.,  101,  4855  (1996). 

29Shafranov,  V.  D.,  Reviews  of  Plasma  Physics,  Vol.  2,  ed.  M.  A.  Leontovich,  p.  103,  Con¬ 
sultants  Bureau,  New  York  (1966). 

30Garren,  D.  A.,  and  Chen,  J.,  Phys.  Plasmas,  1,  3425  (1994). 

3101medo,  O.,  and  Zhang,  J.,  Astrophys.  J.,  718,  433  (2010). 

32Xue,  M.  L.,  and  Chen,  J.,  Solar  Phys.,  84,  119  (1983). 

33Linker,  J.  A.,  Mikic,  Z.,  Lioncllo,  R.,  Riley,  P.,  Amari,  T.,  and  Odstrcil,  D.,  Phys.  Plasmas, 
10,  1971  (2003). 

34Lynch,  B.  J.,Antiochos,  S.  K.,  Li,  Luhman,Y.,  J.  G.,  and  DeVore,  C.  R.,  Astrophys.  J., 
697,  1918  (2009). 

35Karpen,  J.  T.,  Antiochos,  S.  K.,  and  DeVore,  C.  R.,  Astrophys.  J.,  760,  81  (2012). 
36Kadomtsev,  B.  B.,  in  Reviews  of  Plasma  Physics ,  2,  Consultants  Bureau,  New  York,  p. 
178  (1966). 

37Cargill,  P.  J.,  Chen,  J.,  and  Garren,  D.  A.,  Astrophys.  J.,  423,  854  (1994). 

38Torok,  T.,  and  Kliern,  B.,  Astrophys.  J.,  630,  L97  (2005) 


15 


